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ABSTRACT 

We present a description for the angular momentum loss rate due to magnetic braking for late type stars taking into 
account recent observational data on the relationship between stellar activity and rotation. The analysis is based on an 
idealized two component coronal model subject to constraints imposed on the variation of the coronal gas density with 
rotation period inferred from the observed variation of X-ray luminosity, L x , with rotation rate, f2, (L x oc f2 2 ) for single 
rotating dwarfs. An application of the model to high rotation rates leads to a gradual turnover of the X-ray luminosity 
which is similar to the saturation recently observed in rapidly rotating dwarfs. The resulting angular momentum loss 
rate, J, depends on in the form J oc where (3 ~ 3 for slow rotators and ^1.3 for fast rotators. The relation at 
high rotation rates significantly differs from the power law exponent for slowly rotating stars, depressing the angular 
momentum loss rate without necessarily requiring the saturation of the magnetic field. The application of this model 
to the evolution of cataclysmic variable binary systems leads to mass transfer rates that are in approximate accord 
with those observed in comparison to rates based on either a Skumanich law or an empirical law based on [3 — 1. 

Subject headings: binaryrclose - binaryrmagnetic braking - binary:rotation - stars:late-type 



1. INTRODUCTION 

It is generally accepted that the angular momentum 
removed from a stellar surface by the action of a magnet- 
ically coupled stellar wind is important in studies of the 
formation and evolution of close binary systems. Mag- 
netic braking (MB) is thought to be the fundamental 
mechanism responsible for orbital angular momentum 
loss in a number of classes of close binaries, such as 
cataclysmic variables (CVs) and low mass X-ray bina- 
ries (LMXBs). For such systems, MB has two important 
roles. It provides a mechanism for promoting mass trans- 
fer from the low mass donor to its more massive compact 
companion. In addition, it determines the range of post 
common envelope orbital separations for which the sys- 
tem evolves into this mass transfer stage. 

MB was first suggested as a mechanism for remov- 
ing angular momentum from single stars by Schatzman 
(1962), who noticed that slowly rotating stars have con- 
vective envelopes. A key insight was his recognition that 
material lost from the stellar surface is kept in corotation 
with the star by the magnetic field, leading to the result 
that the specific angular momentum carried by the gas 
is significantly greater than in a spherically symmetric 
stellar wind. To maintain corotation, a star should pos- 
sess a substantial magnetic field. This can be achieved 
in low-mass main sequence (MS) stars by gas motions in 
the deep convective envelope, leading to the generation 
of a magnetic field by dynamo action. 

Evidence for angular momentum loss in single stars is 
provided by the observations of young low-mass T Tauri 
stars which have shown that the stellar angular momen- 
tum is significantly higher in the pre-MS stage than in 
the MS stage for stars of similar masses (Stauffer & Hart- 
mann 1987). In addition, the observations of rotational 
velocities of stars in young clusters show a high rate of 
angular momentum loss. Specifically, the K-dwarfs in the 
Hyades (age w 6 x 10 8 yr) have mean rotational veloci- 
ties below 10 km s _1 , but stars of the same spectral type 
in the Pleiades (age « 7 x 10 7 yr) have mean rotational 



velocities of about 40 km s^ 1 (Stauffer 1987). These 
observations supplemented the early work of Skumanich 
(1972) who showed that the equatorial rotation velocities 
of G type MS stars decrease with time, t, as t~ ' 5 . We 
note that this empirical dependence, upon which many 
studies are based, was established for stars with veloci- 
ties in the range of 1 to 30 km s _1 , but its applicability 
to the MS-like companion stars in close binaries of short 
orbital period (where rotational velocities are > 100 km 
s _1 ) is suspect. 

In parallel to these observational developments, many 
theoretical efforts were undertaken to understand the an- 
gular momentum loss mechanism over a wide range in 
rotation rates. For example, the Skumanich relation for 
slowly rotating stars can be reproduced by an angular 
momentum loss rate, J, proportional to f2 3 , where is 
the stellar rotational angular velocity, for a thermally 
driven stellar wind. In this case, the magnetic field is as- 
sumed radial and to vary linearly with f2 (Weber & Davis 
1967). For fast rotation, a wide variety of theoretically 
predicted laws in the form Q oc t~ a are possible with a 
ranging from 0.5 to 4 (Mestel 1984). More recently, Kep- 
pens, MacGregor, & Charbonneau (1995) considered a 
polytropic magnetized solar wind within the framework 
of a model developed by Weber & Davis (1967) show- 
ing that the angular momentum loss rate is consistent 
with the Skumanich law at low angular velocities, but 
the power law dependence is shallower asymptotically 
approaching J oc fl 2 at high angular velocities due to 
centrifugal effects. Additional factors affecting the an- 
gular momentum loss have been pointed out by Mestel 
& Spruit (1987, hereafter MS87) who showed that, in a 
two-component model, more gas is confined in a dead 
zone, which does not contribute to the braking, at high 
angular velocities. 

To further complicate the picture, significant observa- 
tional developments have increased the disparity between 
observations and interpretations based on the Weber- 
Davis type prescription (for a review see Collier Cameron 
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2002). In particular, the observations of stars in open 
clusters have shown that the Skumanich law overesti- 
mates the spin down rate to an age of about 10 8 yr and, 
thus, docs not explain the presence of fast rotators in 
the Pleiades (Stauffer 1987; Andronov et.al 2003). As a 
possible resolution to this problem, MacGregor & Bren- 
ner (1991) suggested that magnetic braking is reduced 
at high rotation rates due to the possible saturation of 
a stellar dynamo resulting in an angular momentum loss 
rate proportional to SI. Here, the strength of magnetic 
field at the stellar surface becomes independent of the 
angular velocity in the convective zone at high rotation 
rates. The saturation regime is normally considered to 
occur at rotational rates of about 10f2©, where f2© is the 
angular velocity of the Sun. While the magnetic field 
strength in active regions of the Sun is a few orders of 
magnitude higher than the average, the fundamental rea- 
son for saturation of the magnetic field at the rotational 
velocity which is only one order of magnitude higher than 
in the Sun is not understood. Another possible mecha- 
nism presumes that the depression of the rate of angular 
momentum loss may be related to the appearance of an 
increasingly complex field topology as the rotation rate 
increases leading to a reduction in the fraction of open 
field lines (e.g., Taam & Spruit 1989). 

Observations of CVs also do not provide direct support 
for an angular momentum loss rate proportional to a 
linear function of the angular velocity (see Andronov ct 
al. 2003). The empirically determined mass transfer 
rates in CVs have been estimated to follow the relation 

logM = 3.31ogP h -H.2±l , (1) 

where M is the mass transfer rate in M© yr -1 and Ph 
is the binary orbital period in hours respectively (Pat- 
terson 1984). Although the inferred mass transfer rates 
are somewhat uncertain, they lie systematically higher 
than one predicts from the application of the angular 
momentum loss rate associated with saturated magnetic 
braking. 

It is therefore clear that there is further need to inves- 
tigate MB for fast rotators. Accordingly, we attempt to 
provide additional insight into the angular momentum 
loss rate associated with the MB process in the context 
of the idealized MS87 model in conjunction with recent 
observations for the X-ray emission of rotating dwarfs 
(Pizzolato et al. 2003). With these observational con- 
straints, we present a modified angular momentum loss 
description and apply it to evolutionary calculations in 
a CV-type close binary system. In §2, the model is de- 
scribed and its application to the evolution of CVs is 
presented. We summarize our results and discuss the 
possible implications of the modified angular momentum 
loss rate in the final section. 

2. MODEL OF MB 

We consider MB within the framework of the two- 
component model described in MS87. In contrast to 
the earlier model described by Weber & Davis (1967), 
a dipole magnetic field, in addition to a radial field, and 
centrifugal acceleration effects on the wind arc included. 
The strong field regions of this dipole field are envisioned 
to confine hot plasma in a corotating dead zone (Mestel 
1968) . Although the X-ray emitting corona is largely pro- 



duced by matter confined in loops (Vaiana, Krieger, & 
Timothy 1973) which fluctuate greatly in the Sun's mag- 
netosphere (Sheeley & Golub 1979; Shimizu & Tsuneta 
1977), we shall assume that the average behavior in this 
time dependent structure can in the lowest order be ap- 
proximated by a steady hydrostatic dead zone model. 
Partial support for such a simple two-component coronal 
structure is provided by the observations taken during so- 
lar eclipse using the Skylab telescope and Yohkoh X-ray 
satellite. Since the matter in the dead zone is trapped 
within the magnetosphere and is not lost in the wind, 
the efficiency of MB in this idealized model is reduced in 
comparison to other idealized theoretical models without 
such regions. The dipole field in the model is assumed 
to be closed up to a cusp point that defines the radius 
of the dead zone, i?d- Beyond this point the magnetic 
field is assumed to be essentially radial. Recently, sup- 
port for the large scale description of such a picture has 
been provided by the X-ray observations of the Sun ob- 
tained from the Ulysses mission. In particular, Li (1999) 
showed that the assumption of homogeneity in latitude 
of the wind flux and the radial magnetic field at large dis- 
tances, a basic feature inherent in the model, is basically 
correct. 

2.1. Constraints on the coronal density from L x 

It was noted in MS87 that the two-component mag- 
netic field model should be consistent with the X-ray 
observations of single stars. Recent results by Pizzolato 
et al. (2003) confirm earlier results by Pallavicini et al. 
(1981), showing that the X-ray luminosity, L x , is propor- 
tional to SI 2 for slowly rotating stars. Since the X-ray 
emission is related to the coronal gas density in the dead 
zone, po.d, in this model, /9o,d should be a function of SI. 
In the simplest model consistent with the X-ray observa- 
tions, MS87 considered the possibility that po,d is linear 
in SI. This choice was based on the assumption that the 
observed variation of coronal emission is determined pri- 
marily by the variation in density, that is, L x oc p\ d and 

, hence, L x oc SI 2 , neglecting the weak dependence on 
coronal temperature. Such a picture is an oversimplifi- 
cation since the variation in X-ray luminosity not only 
reflects variations in the coronal gas density, but also the 
volume of the X-ray emitting zone. 

Further compounding the validity of the linear relation 
between the coronal gas density and rotational angular 
velocity are the recent observations of X-ray emission 
from late-type dwarfs, which have shown that the rela- 
tionship between the X-ray luminosity and rotation rate 
changes form at high rates (Pizzolato et al. 2003). In 
particular, the X-ray luminosity is found to be nearly 
independent of SI for SI > Sl x ~ 2 — 12f2© in the mass 
range of 0.22 — 1.29 M©. The qualitative change from 
L x oc SI 2 to L x independent of SI (for O > Q x ) provides 
an important constraint on the MB model. In the fol- 
lowing we make use of this observed trend to constrain 
the variation of po,d with SI in the MS87 model. 

In order to estimate the X-ray luminosity from the 
MS87 model, we assume that the emitting region can 
be roughly identified with the dead zone in some spa- 
tially and time averaged sense. The radius of the dead 
zone is taken from eq. (8)ms i n MS87, where it is as- 
sumed that p ,d oc S2 P and B oc SI. In addition, we 
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Fig. 1. — The variation of the X-ray luminosity and angular 
momentum loss rate for a 1 Mq star as a function of angular ve- 
locity. Upper panel: The variation of L x , as a function of CI. L x .q 
corresponds to L x at CI = CIq. The thick solid line represents the 
solution for which p oc CI 06 . Thick dash-dotted line represent the 
solution for the same density variation but with the saturation of 
the magnetic field at CI > 10^0. The upper and lower thin dash- 
dotted lines show L x for p oc Cl° and p oc CI respectively. Lower 
panel: The variation of J. as a function of CI Jq corresponds to 
the rate of the angular momentum loss for the star rotating with 
CI = CIq. Lines as described for the upper panel. 



assume that the parameter £d = BQ/Sirpo^a^ is equal 
to 60 where ad is the sound speed in the dead zone (see 
MS87). For a solar type star of mass, M, equal to 1 M 
the parameter p was varied such that the L x varied with 
Q, in approximate accordance to the observed relation- 
ship, with a normalization chosen to correspond to the 
solar X-ray luminosity at CIq. The functional form of 
L K with respect to fi was found to be very sensitive to 
p for f2 < IOOq, varying as fi 1-2 and f2 2 8 for p equal to 
and 1 respectively. The relation between L x and f2 is 
illustrated in the upper panel of Fig. 1, where the best 
fit for L x oc il 2 corresponds top« 0.6. 

Fig. 1 also reveals that L x is relatively insensitive to 
angular velocity at high angular velocities, especially for 
p > 0.6. This tendency toward saturation occurs at 
high rotation rates since the centrifugal effects become 
increasingly important, resulting in the shrinkage of the 
dead zone (MS87). The decrease in volume tends to off- 
set the increase in p ad as Q is increased resulting in a 
nearly constant value for L x . For this model, i x actu- 
ally decreases for sufficiently large Q. At larger values 
of p, L x increases monotonically since the density de- 
pendence more than compensates for the reduction in 
emitting volume. On the other hand, the volume effect 
is more important for smaller values of p, leading to a 
significant reduction in i x as fi is increased. 

The ratio of i x at in the saturated regime, taken to 
correspond to the maximum value of L x , to L XjQ at £1 = 
Cl & for p — 0.6 is in approximate accord to that observed. 
Specifically, A = logi x at /i x . Q = 2.3 from the numerical 
results ior p = 0.6 and A — 2.4 from the observations (see 
Fig. 5 in Pizzolato et al. 2003). We point out that the 



qualitative behavior of L x as a function of il is insensitive 
to the adopted choice of Quantitatively, the detailed 
results only change slightly. For example, for Q = 4 (see 
MS87), p ~ 0.4 and A = 2.2. 

For comparison, we have also displayed the solution 
with a saturated magnetic field for f2 > 10 CIq in Fig. 
1. Although L x is seen to turnover as well, we note that 
it occurs at a level significantly less than observed with 
A = 2. We cannot exclude the possibility that saturation 
of the magnetic field takes place based on the observed 
relation L x (Cl) since it may take place at rotation rates 
higher than IOCIq making it difficult to distinguish the 
influence of the saturation of the magnetic field from the 
influence of the variation of the dead zone with angular 
velocity. However, the X-ray observations of solar type 
stars alone do not require the additional assumption of 
a saturated magnetic field at high rotation rates. 

2.2. MB rate 

The transport of the angular momentum by the wind 
and by Maxwell stresses is equivalent to that found by 
assuming corotation maintained out to the Alfven surface 
(Mcstel 1984). The angular momentum loss rate can be 
written as 

,77/2 

-J = 4ttO / (p A v A R 2 A ){R A sin 6) 2 sin 9d9, (2) 
Jo 

where 9 is the polar angle, and p A , v A , and R A are the 
density, velocity, and radius at the Alfven surface respec- 
tively. By continuity, p A v A /p d v d = B A /B d = (R d /R A ) 2 
assuming a radial field beyond the dead zone region. 
Here the radius od the dead zone is found from equ. (8)ms 
of MS87. Then 

-<*&)(£)''("■£)■ < 3 > 

where R is the stellar radius and Cj incorporates all 
physical constants. 1 The Alfven radius is found from 
equ. (14) MS of MS87. 

In the lower panel of Fig. 1, the angular momentum 
loss rate relative to the Sun, log(J/j0), is shown as a 
function of normalized angular velocity, CI/CIq, for three 
different values of p corresponding to p = 0, 0.6 and 1. 
For reference, we also illustrate the variation of the an- 
gular momentum loss rate for the case of a saturated 
magnetic field (for Q > Cl x ) with p = 0.6. For slow rota- 
tion, the MB is consistent with the Skumanich law and 

1 F(M, CI/CIq) also depends on the solar sonic velocity and the 
coronal density, but we treat them as constants in order to scale 
the equation to the observed value of the solar Alfven radius and 
the dead zone. We are aware that R A q and q are not known 
precisely, and in particular are different in MS87 and Li(1999). 
Since we do not attempt to derive the precise value for the current 
solar J directly from observations as in Li (1999), we adopted these 
values as in MS87 and assume that uncertainties in the measured 
solar values arc absorbed in the calibration of J. 
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is nearly independent of the dependence of the coronal 

density on angular velocity (in contrast to the variation Table 1 - The rotational velocities obtained from the theoretical 

c T n\ TJ c c i- 4. j.- /o ^ mr> \ MODELS AND FROM THE OBSERVATIONS FOR THE PLEIADES AND HYADES CLUS- 

of L x with fi). However for fast rotation {Si > 10 « Q ), TERSi IN KM s -i. H ere, "observed" corresponds to the maximum ob- 
J oc SI 1 ' 3 for p = 0.6, the value inferred from the con-SERVED velocity of stars in the Pleiades or the Hyades, "saturated" 
straints based on the X-ray observations (see above). We represents the results of calculations based upon the saturated mag- 

, , , , , ,. f f , , ' NETIC BRAKING LAW, AND "STANDARD" REPRESENTS THE MAXIMUM ROTATIONAL 

note tnat tne asymptotic torm ot tne angular momentum yELOCITY THAT CAN BE OBXAINED WIXH THE STANDARD prescription for the 
loss rate on the angular velocity, in contrast to the X-ray magnetic braking (all the data is taken from Andronov et al. 2003); 
luminosity, is insensitive to the value of p if p lies be- "this work" represents results of the spin-down using equ|3] 
tween 0.6 and 1. On the other hand, if the magnetic 
field is saturated, the angular momentum loss rate is ap- 
proximately constant at high rotation rates. Hence, the 
loss rate from a centrifugally driven wind differs from the 
linear functional form suggested by MacGregor & Bren- 
ner (1991) for saturation of the stellar dynamo based on 
a thermally driven stellar wind model. 

Based on these results, the solution for the rate of an- 
gular momentum loss due to MB can be presented in the 
parameterized form 
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Here, the value of Kj depends on the mass of the star, 
and the temperature and density in the corona which 
can be derived directly from the observational data (Li 
1999). This parametrized form is presented only to pro- 
vide a more transparent form of the angular momentum 
loss rate for comparison between the magnetic braking 
law based on equ. and the power law representations 
describing the Skumanich law (e.g., as in Verbunt & 
Zwaan 1981) or the saturated braking law (Andronov 
et al. 2003). In detailed calculations, the rate of an- 
gular momentum loss due to MB should be found us- 
ing equ. |3| together with highly nonlinear equ. (8)ms an d 
equ. (14)ms of MS87. 

In this study we adopt the approach of calibrating Kj 
using evolutionary calculations to reproduce the known 
solar rotation period at the age of the Sun (see Kawaler 
1988). For Sl Q = 3 x 10~ 6 we find Kj = 6 x 10 30 [dyn cm], 
which is about three times the value obtained in Li (1999) 
and about the same as in Pylyser & Savonije (1988) (their 
value includes a variable factor 0.73 < / < 1.78 ). The 
corresponding value for the constant in the equ. (4) is 
Cj = 2.1 x 10 27 [dyn cm]. 2 

Our description of magnetic braking has been applied 
to a study of the rotational evolution of single stars for 
comparison to the observed spin-down of stars in young 
clusters. Although stellar rotation is not included in the 
stellar structure, we calculated two sets of stellar models 
of mass 1, 0.8 and 0.6 M Q with the chemical composi- 
tion as in Pleiades and Hyades to obtain an indication 
of the effect of our modified braking rate. The models 
were evolved for a time corresponding to 7 x 10 7 yr for 
the Pleiades and 6 x 10 8 yr for the Hyades clusters. The 
initial surface rotational velocity was taken as the break- 
up velocity. The numerical results as presented in Table 

2 Our code does not incorporate the angular velocity in the equa- 
tion of hydrostatic equilibrium and therefore a more accurate cali- 
bration can be carried out with an evolutionary code including the 
effects of rotation. However, since the value of Kj mainly depends 
on the evolution of a relatively slowly rotating star, where Q does 
not contribute significantly to hydrostatic equilibrium, our results 
should be adequate. 



1 show that the level of agreement between the observa- 
tions and the angular momentum loss rate presented in 
this study and with the loss rate based on a saturated 
model is comparable. The agreement is good for the 
Pleiades cluster and is slightly better for stars of mass 
1 and 0.8 M in comparison to the lower mass model 
for the Hyades cluster. Here, the observational data for 
v sin i and data for the spin-down based on the saturated 
magnetic braking law and the Skumanich law have been 
taken from Andronov et al. (2003). 

2.3. Binary evolution 

To determine the consequences of our modified angu- 
lar momentum loss rate prescription on the evolution of a 
binary system, we calculated the mass transfer phases of 
two binary model sequences. The model stars are chosen 
to be of a solar metallicity. The first sequence consisted 
of 1 Mq evolved donor and the second a 0.8 Mq un- 
evolved donor, both are on the MS and, with a 0.6 M & 
white dwarf companion. Such types of systems can be 
considered to be representative of CVs. For the evolved 
donor, the central hydrogen abundance was X = 0.35. 
The radii of the evolved and unevolved model stars were 
such that they filled their Roche lobes at orbital periods 
of P = 7.5 h and 5.2 h respectively. To follow the evolu- 
tion, we used the code described in Ivanova et al. (2003), 
supplemented with angular momentum loss rates associ- 
ated with MB based on equ. (3). The evolution of the 
system is assumed to be non conservative with nearly all 
the mass lost from the donor 3 ejected as a consequence of 
nova explosions or as a result of an optically thick wind 
from the white dwarf surface (Hachisu, Kato, & Nomoto 
1999). For comparison, we also performed mass trans- 
fer sequences with MB according to Verbunt & Zwaan 
(1981) based on the Skumanich law (taken from Pylyser 
& Savonije, 1988). 

The results for the mass transfer rates are shown as a 
function of orbital period in Fig. 2. The average mass 
transfer rate for both donors is ~ 5 x 10 -10 M Q yr -1 . For 
comparison, we also display the mass transfer rates cor- 

3 For details on the treatment in our code for the WD evolu- 
tion and its ability to accumulate mass see Ivanova & Taam in 
preparation (see also Li & van den Heuvel 1997). 



5 



-8 

u 

>* 

© 

S 

OB 

o 

~ -9 



1 1 

1 f 

o.8M s , J«n3 1 f I 

! ! 
\ i 

\ / 


\ 

\ 

\ 

\ 


/ 1 

im , | / \ 

! 






0.8M S , j«n L3 







I & I I I I I I I I I I I 

2 4 6 8 

P[h] 



Fig. 2. — The variation of mass transfer rate (in Mq yr — 1 ) as 
a function of orbital period (in hours) for two model sequences. 
The evolution of a 1 Mq evolved donor and 0.6 Mq white dwarf is 
illustrated by a dotted line for the angular momentum prescription 
presented in this paper and as a dash dashed line for the angular 
momentum prescription assuming that the Skumanich law applies 
to high angular velocities. The dot dashed and dashed line il- 
lustrate the evolution of a 0.8 Mq unevolvcd donor with a 0.6 Mq 
white dwarf with the prescription adopted here and the Skumanich 
law respectively. For reference, the solid line represents the empir- 
ically inferred mass transfer rates obtained by Patterson (1984). 



responding to those observationally estimated according 
to an empirical relation obtained by Patterson (1984). 
It can be seen that the mass transfer rates promoted by 
the Skumanich law are generally higher than the empir- 
ical values, and the rates produced by the modified MB 
rate are in approximate agreement. The evolutionary 
timescale of the systems is increased by about a factor of 
5 for the both sequences with the modified MB rate in 
comparison to the prescription based on the Skumanich 
law. We note that this is significantly less than the factor 
of 100 found by Andronov et al. (2003) in their case for 
MB rate based on J oc ft. 

3. SUMMARY 

The X-ray luminosity of late type stars and the angular 
momentum loss rate associated with magnetic braking of 
these stars has been examined in terms of an idealized 
two-component coronal model in the light of recent ob- 
servational studies relating stellar activity to stellar ro- 
tation. By assuming that the X-ray luminosity can be 
used as a proxy for the coronal gas density in the mag- 
netosphere we have, as a first step, constrained the two- 
component model to reproduce the observed variation 
of X-ray luminosity with stellar rotation (L x oc f2 2 ) for 
slowly rotating late type dwarfs. The application of this 
model to rapidly rotating dwarfs provides a qualitative 
explanation for the observed saturation at high rotation 
rates where the decrease in volume of the X-ray emit- 
ting coronal region nearly compensates for the increase 
in coronal gas density with increasing rotational velocity. 

The angular momentum loss rate resulting from such 
a model has been analyzed and a parameterized form for 
the rate has been presented. In particular, it has the form 



J oc fi 3 for slow rotators to reproduce the Skumanich law 
and J oc f2 13 for fast rotators (f2 > 10f2©). The angular 
momentum loss rate is more sensitive to SI in the high 
rotation rate regime than the form, J oc 17, adopted in 
recent studies by Andronov et al. (2003). In addition, 
the reduction of the angular momentum loss rate from 
that based on the Skumanich relation for slowly rotat- 
ing stars, required by the observation of rapidly rotating 
stars in clusters, can be accommodated without neces- 
sarily invoking the assumption that the magnetic field 
saturates at high rotation rates. We note, however, that 
the observations of rotational vclocites of stars less mas- 
sive than about 0.4 M Q (see Andronov et al. 2003 and 
references therein) suggest that the angular momentum 
loss rate is further weakened. 

The modified form of the angular momentum loss rates 
can be applied to both LMXBs and CVs, and we have 
carried out binary evolutionary calculations for model 
systems representative of CVs as an illustration. It has 
been shown that the mass transfer rates are approxi- 
mately consistent with those observed for CVs provided 
that the observed rates are accurate and reflect the long 
term secular values. The calculated rates for unevolved 
and evolved main sequence-like stars are found to be 
about an order of magnitude lower than rates obtained 
using angular momentum loss rates based on the appli- 
cation of the prescription for slowly rotating stars to the 
rapidly rotating stars in CVs and are higher by a fac- 
tor of 10 than the rates obtained using rates obtained 
using a saturated braking law (Andronov et al. 2003). 
As a consequence the time scale for both the pre CVs 
to evolve into the CV phase and the evolution during 
the CV phase are increased in comparison to previous 
detailed binary evolutionary studies, but a factor of 10 
shorter than obtained using a saturated law. Hence, the 
possibility that finite age effects are important in under- 
standing the distribution of such systems (e.g., King & 
Schenker 2002; Andronov et al. 2003) is made less likely. 

Although the mass transfer rates are significantly 
higher than that obtained by Andronov et al. (2003), 
the rates are a factor of 2 to 3 times lower than the 
rates (M > 1.5 — 2 x 10~ 9 M yr _1 ) required to force an 
unevolved donor sufficiently out of thermal equilibrium 
to explain the period gap width in the disrupted mag- 
netic braking model (see McDermott & Taam 1989; Kolb 

2002) . As noted previously, the upper boundary of the 
period gap at about 3 hours approximately corresponds 
to a main sequence star of ~ 0.4 M©. However, there 
is no observational evidence indicating a discontinuous 
change in the rate of spin down of single stars, although, 
for less massive stars, the rotational velocity data of clus- 
ters suggests a further weakening of the angular momen- 
tum loss rate (see Sills, Pinsonneault, & Terndrup 2000). 
This suggests that other mechanisms should be sought 
to prevent systems from entering into the period gap. A 
possible alternative involves a period bounce first sug- 
gested by Eggleton (1983) involving additional angular 
momentum losses, perhaps, beyond that considered in 
the present study (see e.g., Taam, Sandquist, & Dubus 

2003) . 

An additional consequence of the reduced angular mo- 
mentum loss rate for rapidly rotating stars will be a re- 
vision of the binary orbital period delineating those sys- 
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terns which evolve to shorter periods (angular momen- 
tum loss dominated evolution) from those which evolve 
to longer orbital periods (nuclear evolution dominated). 
This period not only depends on the rate of angular mo- 
mentum loss, but is also a function of the mass of the 
binary components and the degree to which mass is lost 
from the system. For the case of conservative mass trans- 
fer, Pylyser & Savonije (1988, 1989) estimate a bifurca- 
tion period of about 12 hours for systems with an accret- 
ing compact companion. With a reduced effectiveness of 
magnetic braking, it is expected that the transition will 
shift to shorter orbital periods. 

Future phcnomcnological studies of the X-ray emission 
from late type stars and the angular momentum losses 
associated with magnetic braking in these stars are de- 
sirable to determine the explicit mass dependence on the 



input parameters (magnetic field, density and tempera- 
ture of the magnetically confined region) in the simpli- 
fied magnetically coupled stellar wind models. The con- 
straints imposed by both the X-ray studies and rotational 
velocity studies may provide insight into the further re- 
duction of the angular momentum loss rates for stars less 
massive than 0.4 M . The incorporation of these ideas in 
binary evolutionary sequences will be necessary in order 
to make meaningful comparisons between the observed 
and theoretically predicted CV period distributions in 
future population synthesis studies. 
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